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DECOUPLED MAGNETOQUASISTATIC 
NON-LINE-OF-SIGHT POSITION AND 
ORIENTATION SENSING FOR ARBITRARY 
DISTANCES 


CROSS REFERENCE TO RELATED 
APPLICATIONS 


The present application claims priority to U.S. Provisional 
Patent Application No. 62/060,714, filed on Oct. 7, 2014, US 
Provisional Patent Application No. 62/079,738, filed on 
Nov. 14, 2014, and U.S. Provisional Patent Application No. 
62/160,129, filed on May 12, 2015, and is related to U.S. 
patent application Ser. No. 12/824,894, filed on Jun. 28, 
2010, the disclosure of each of which is incorporated herein 
by reference in its entirety. 


STATEMENT OF INTEREST 


The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 
Contractor has elected to retain title. 


TECHNICAL FIELD 


The present disclosure relates to non-line-of-sight posi- 
tioning. More particularly, it relates to magnetoquasistatic 
non-line-of-sight positioning. 


BRIEF DESCRIPTION OF DRAWINGS 


The accompanying drawings, which are incorporated into 
and constitute a part of this specification, illustrate one or 
more embodiments of the present disclosure and, together 
with the description of example embodiments, serve to 
explain the principles and implementations of the disclosure. 

FIG. 1 illustrates exemplary configurations of the two- 
axis transmit (105) and receive (110) systems. 

FIG. 2 illustrates a 2-axis transmitter and receiver system 
during range and orientation experiments on a building 
rooftop. 

FIG. 3 illustrates measurements of the orientation-invari- 
ant range function compared to theory. 

FIG. 4 illustrates the error in range and orientation solu- 
tion. 

FIG. 5 illustrates measurements of the range-invariant 
orientation function compared to theory. 

FIG. 6 illustrates the orientation error and range inversion 
for measurements with variation in orientation at a fixed 
range. 

FIG. 7 illustrates transmit-receive coil coupling. 

FIG. 8 illustrates a two-axis orthogonal transmit and 
receive system. 

FIG. 9 illustrates a transmit and receive system for 
experimental setups. 

FIGS. 10-13 illustrate measurements of the azimuth- 
direction function compared to theory for variation in azi- 
muth-direction. 

FIG. 14 illustrates a two-dimensional map of measure- 
ments at ranges of 10 m, 20 m, and 30 m, and for various 
transmit orientation measurements at a range of 20 m. 

FIG. 15 illustrates a distance experiment setup and mea- 
surements. 

FIG. 16 illustrates the estimated distance error using 
complex image theory, and for the exact integral solutions. 
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FIG. 17 illustrates the 3D problem with 3-axis transmit- 
ting coils and their respective magnetic moments. 

FIG. 18 illustrates a mathematical description of the 3D 
position sensing problem. 

FIG. 19 shows the elevation-direction angle ratio func- 
tion. 


SUMMARY 


In a first aspect of the disclosure, a method is described, 
the method comprising: providing a receiver, the receiver 
comprising a first and second receiver coils configured to 
receive quasistatic magnetic fields in a first and second 
mutually orthogonal directions transmitted by a transmitter, 
the transmitter comprising a first and second transmitter 
coils and configured to transmit the quasistatic magnetic 
fields in a third and fourth mutually orthogonal directions; 
receiving the quasistatic magnetic fields through the 
receiver; measuring, with the receiver, a first voltage origi- 
nating from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the first receiver coil; 
measuring, with the receiver, a second voltage originating 
from the quasistatic magnetic fields due to coupling between 
the second transmitter coil and the first receiver coil; mea- 
suring, with the receiver, a third voltage originating from the 
quasistatic magnetic fields due to coupling between the first 
transmitter coil and the second receiver coil; measuring, 
with the receiver, a fourth voltage originating from the 
quasistatic magnetic fields due to coupling between the 
second transmitter coil and the second receiver coil; based 
on the first, second, third, and fourth voltages, calculating a 
distance between the receiver and the transmitter, wherein 
the measuring of each of the first, second, third, and fourth 
voltages is independent from the remaining voltages. 

In a second aspect of the disclosure, a method is 
described, the method comprising: providing a receiver, the 
receiver comprising a first and second receiver coils con- 
figured to receive quasistatic magnetic fields in a first and 
second mutually orthogonal directions; providing a trans- 
mitter, the transmitter comprising a first and second trans- 
mitter coils and configured to transmit the quasistatic mag- 
netic fields in a third and fourth mutually orthogonal 
directions; transmitting the quasistatic magnetic fields 
through the transmitter; receiving the quasistatic magnetic 
fields through the receiver; measuring, with the receiver, a 
first voltage originating from the quasistatic magnetic fields 
due to coupling between the first transmitter coil and the first 
receiver coil; measuring, with the receiver, a second voltage 
originating from the quasistatic magnetic fields due to 
coupling between the first transmitter coil and the second 
receiver coil; measuring, with the receiver, a third voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the first 
receiver coil; measuring, with the receiver, a fourth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the second 
receiver coil; based on the first, second, third, and fourth 
voltages, calculating a distance between the receiver and the 
transmitter, wherein the measuring of each of the first, 
second, third, and fourth voltages is independent from the 
remaining voltages. 

In a third aspect of the disclosure, a method is described, 
the method comprising: providing a receiver, the receiver 
comprising a first, second and third receiver coils configured 
to receive quasistatic magnetic fields in a first, second and 
third mutually orthogonal directions transmitted by a trans- 
mitter, the transmitter comprising a first, second and third 
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transmitter coils and configured to transmit the quasistatic 
magnetic fields in a fourth, fifth and sixth mutually orthogo- 
nal directions; receiving the quasistatic magnetic fields 
through the receiver; measuring, with the receiver, a first 
voltage originating from the quasistatic magnetic fields due 
to coupling between the first transmitter coil and the first 
receiver coil; measuring, with the receiver, a second voltage 
originating from the quasistatic magnetic fields due to 
coupling between the first transmitter coil and the second 
receiver coil; measuring, with the receiver, a third voltage 
originating from the quasistatic magnetic fields due to 
coupling between the first transmitter coil and the third 
receiver coil; measuring, with the receiver, a fourth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the first 
receiver coil; measuring, with the receiver, a fifth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the second 
receiver coil; measuring, with the receiver, a sixth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the third 
receiver coil; measuring, with the receiver, a seventh voltage 
originating from the quasistatic magnetic fields due to 
coupling between the third transmitter coil and the first 
receiver coil; measuring, with the receiver, an eighth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the third transmitter coil and the second 
receiver coil; measuring, with the receiver, a ninth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the third transmitter coil and the third 
receiver coil; based on the first, second, third, fourth, fifth, 
sixth, seventh, eighth and ninth voltages, calculating a 
distance between the receiver and the transmitter, wherein 
the measuring of each of the first, second, third, fourth, fifth, 
sixth, seventh, eighth and ninth voltages is independent from 
the remaining voltages. 

In a fourth aspect of the disclosure, a method is described, 
the method comprising: providing a receiver, the receiver 
comprising a first, second and third receiver coils configured 
to receive quasistatic magnetic fields in a first, second and 
third mutually orthogonal directions; providing a transmit- 
ter, the transmitter comprising a first, second and third 
transmitter coils and configured to transmit the quasistatic 
magnetic fields in a fourth, fifth and sixth mutually orthogo- 
nal directions; receiving the quasistatic magnetic fields 
through the receiver; measuring, with the receiver, a first 
voltage originating from the quasistatic magnetic fields due 
to coupling between the first transmitter coil and the first 
receiver coil; measuring, with the receiver, a second voltage 
originating from the quasistatic magnetic fields due to 
coupling between the first transmitter coil and the second 
receiver coil; measuring, with the receiver, a third voltage 
originating from the quasistatic magnetic fields due to 
coupling between the first transmitter coil and the third 
receiver coil; measuring, with the receiver, a fourth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the first 
receiver coil; measuring, with the receiver, a fifth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the second 
receiver coil; measuring, with the receiver, a sixth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the second transmitter coil and the third 
receiver coil; measuring, with the receiver, a seventh voltage 
originating from the quasistatic magnetic fields due to 
coupling between the third transmitter coil and the first 
receiver coil; measuring, with the receiver, an eighth voltage 
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originating from the quasistatic magnetic fields due to 
coupling between the third transmitter coil and the second 
receiver coil; measuring, with the receiver, a ninth voltage 
originating from the quasistatic magnetic fields due to 
coupling between the third transmitter coil and the third 
receiver coil; based on the first, second, third, fourth, fifth, 
sixth, seventh, eighth and ninth voltages, calculating a 
distance between the receiver and the transmitter, wherein 
the measuring of each of the first, second, third, fourth, fifth, 
sixth, seventh, eighth and ninth voltages is independent from 
the remaining voltages. 


DETAILED DESCRIPTION 


Existing radiolocation systems such as the Global Posi- 
tioning System (GPS), Radio-Frequency Identification 
(RFID) systems, and Ultra Wide-Band (UWB) systems are 
based on propagating electromagnetic (EM) waves. These 
systems have reduced accuracy in non-line-of-sight (NLoS) 
environments due to propagation losses, delays, or multi- 
path effects, see for example Refs. [4, 6, 7]. In contrast, 
magnetoquasistatic (MQS) fields are not significantly dis- 
turbed and can offer accurate radiolocation in NLoS envi- 
ronments, see Ref. [3]. Magnetoquasistatic fields are a class 
of electromagnetic fields in which a slowly oscillating 
magnetic field is dominant. A magnetoquasistatic field is 
typically generated by low-frequency induction from a mag- 
netic dipole or a current loop or coil. In the past, MQS 
systems using frequencies of up to a few kiloHertz were 
developed for above ground positioning, see Ref. [8, 9]. 
These techniques did not account for induced eddy-currents 
in the ground, which limited the range of operation to a few 
meters above the ground. However, this simplification in the 
technology did permit a decoupled solution for position and 
orientation, see Ref. [9]. The decoupling of position and 
orientation sensing is important in magnetoquasistatic sys- 
tems to provide linear unambiguous solutions and also for 
low-latency or fast solution times. 

The solutions are decoupled when both the position and 
orientation of an object can be determined through equations 
that are not coupled and can be solved separately and 
independently. Long-range MQS positioning above ground 
was recently demonstrated using frequencies of a few- 
hundred kiloHertz and by accounting for the induced eddy- 
currents in the ground, see Ref [3]. However, the strong 
ground effects did not permit a decoupled solution for range 
and orientation of the device to be radiolocated, and instead 
required solution of non-linear coupling equations. Due to 
use of the non-linear field coupling equations, see Ref. [3], 
previous attempts to extend the long-range MQS technique 
to 2D and 3D have resulted in solutions that exhibit large 
errors (1-3 m), which are strongly varying as a function of 
the orientation and position of the transmitter and receivers, 
see Ref. [2]. Furthermore, the complexity in the non-linear 
solution space impedes the ability to converge to a correct 
solution, or to solve in real-time. Attempts to use orthogonal 
fields to aid the non-linear convergence has resulted in slight 
improvements in accuracies at the expense of significantly 
more complicated algorithms with high computational com- 
plexity, see Ref [1]. 

The present disclosure describes a long-range MQS tech- 
nique to accurately determine the orientation-invariant range 
and range-invariant orientation of a device above ground. To 
achieve this, the non-linear field equations describing the 
coupling, which includes ground effects, are decoupled into 
range and orientation field equations that can be directly 
inverted. 
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FIG. 1 illustrates exemplary configurations of the two- 
axis transmit (105) and receive (110) systems. The trans- 
mitted fields generate eddy-currents which create secondary 
fields. These secondary fields are approximately modeled 
using complex image theory for simplicity, see Ref. [3]. FIG. 
1 illustrates also a top view of the systems (115). 

The vector magnetic field of an electrically-small coil 
driven by a current is given by (see Ref. [3]): 


(1) 


Bey, y, Z)= Ho =}, 


An P 


where r =r? is the range vector centered at the origin of 
the coil, m is the vector magnetic moment of the transmitter 


and where the moment is assumed unity, i.e., Iml=1 Am?. 
The voltage measured by a receiving coil within the quasi- 


static region is given by V=-jo[B fla, where a is the surface 
area of the receive coil, and 11 is the receiver’s surface 
normal unit vector. Induced eddy-currents in the ground 
generate secondary fields detected by the receiver, which 
can be modeled using complex image theory as indicated in 
FIG. 1 (see Ref. [3]). The total voltage measured at the 
receiver is given by V7=V,,,_.,,+Vini+n, Where V,,_.,, is the 
voltage induced in the receiver due to the transmit direct 
fields, V,,,;_.,, 18 the voltage induced in the receiver due to the 


image fields of m’. The total voltage is therefore given by: 


(2) 


jblod| 3ih-PA-P)— th-f | 3G Pi MA) — fh 


Vr = 
Te Say rn a 


Of concern are the radial ranging near the 2D xy-plane, 
and the orientation of the transmitter about the Z axis. 
Horizontal magnetic moments are chosen for the transmit 
and receive (vertical coils), and cylindrical coordinates are 
used to describe the problem. In the cylindrical system, 
f=(p/r)p+(z /r)Z, where r=(p?+z7)'”?, and #.=(p/r,)p+[(z+H)y/ 
r,]Z, where r=(p?+(z+H)?)'”, and H=2h+a.. Here a=8(1-j), 
where 6 is the skin depth of penetration in the ground. Due 
to the use of horizontal magnetic moments or vertical coils 
(that implies 0=90°), -Z=1h,-z=f-z=,-7=0, and j=n=m,, 
which substantially simplifies the total voltage due to the 
transmitter and ground to: 


(3) 


Oya 3(plry (hh PM-) ah im 3(elr;y (PAP) —My-h 
An n : 


i 


pr 


Due to operation close to the 2D xy-plane, it is assumed 
that z—0 and hence p/r=1. The asymptotic convergence 
from the source field at near ranges, to the complex image 
contribution at far ranges, is understood in the context of 
complex image theory, see Ref [3]. At short ranges (o<<H), 
it is found that p/r~p/H-0 and 1/r?>>1/r7, so that the 
source field significantly dominates the sensed voltage. On 
the other hand, at long ranges (p>>H), it is found that 
p/r,1 and 1/r~1/r,7, so that both source and image fields 
contribute to the sensed voltage, but the scaling of the 
dot-product in the numerator of the image contributions are 
not significant. In summary, at short ranges only the source 
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field is significant, but p/r~1, whereas at long ranges, both 
source and image fields are important, but p/r,1. The total 
voltage can now be approximately re-written as: 


Oflya (4) 


An 


Vp x —j 2 [30h -p)-p) —h-f] x] > + = 
r i 


where the r and r, are kept in the denominator since the 
cubic dependence offers high sensitivity to z or H. The term 
on the right (1/r°+1/r,3), is the geometric distance function 
from complex image theory, see Ref. [3], and can be 
replaced by using higher-order discrete complex images (see 
Ref. [5]) when slightly higher accuracy is needed at pH. 
However, the higher-order images provide diminished 
returns due to further depths. 

The description of the voltage sensed at the receiver given 
in Eq. (4) is simpler than that used in previous documents, 
see Refs. [2, 3]. However, it is still non-linear, as the 
orientation of the transmitter, orientation of the receiver, and 
range functions are strongly coupled. Previous attempts to 
extend Eq. (2) to 2D and 3D (see Ref [2]) have resulted in 
large errors (1-3 m), which are also strongly dependent on 
the orientation and position of the devices. The non-linear 
solution space is also complex and often prevents conver- 
gence to a correct solution or any real-time solution. 

To solve this problem, it can be noted that Eq. (4) has a 
similar form to the traditional free-space magnetoquasistatic 
(MQS) coupling equation, and thus can be decoupled in a 
related manner to that in short-range MQS techniques, see 
Ref. [9]. A 2-axis transmit and receiver concept as shown in 
FIG. 1 is used, where each axis (1h,, ,) of the transmitter 
system, and similarly each axis (f,, f,) of the receiver 
system, are horizontal and orthogonal to each other. The 
orientations of the transmitter and receiver systems are 
arbitrary with respect to each other. Due to this 2x2 con- 
figuration, there is a total of four field equations generated. 
These four equations can be considered by first focusing on 
the contributions at both receivers due to transmitter 1. At 
receiver 1, taking note that p-1h,=cos ,,,, P fi,=Cos $,,, 
and 1h,-f,=cos(q,,;-$,,, )=Cos 9,,;+8iN >,, sin ,,;, it is 
possible to re-write the total voltage for this case as: 


1 1 
BT Rl 


For receiver 2, similar steps can be carried, while noting 
that 9,.=$,,,+7/2, which allows a re-write of the total 
voltage: 


(5) 


Wfloa 
Vromtont) = —F An 


[2cosn1 COSPrnt — SIND 1 SIND m1 )] X 


Woda 
eT 


Vromion2) = —J 7 —2cos$ni SING) — COSHn1 SIND nt )] X | | 


Now upon taking the squares of the voltage descriptions 
and summing both receiver contributions, it can be found 
that: 


(7) 


2 
VF¢ntnl) 


-{% 
+ Vironton2) 


where c=-j@u,o/4a is a constant. A point to be noted here 
is that regardless of the orientation of m,, the orientation of 
the receivers do not have an effect in Eq. (7). By taking 


= eo + 3cos(2bm1) E . |: 
2 PP G 
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similar steps, the counterpart to Eq. (7) that describes the 
case of the transmitter m, can be found: 


ha =| 


where the change 9,,.=$,,;+7/2 was made so that both 
Eqs. (7) and (8) have the same variables. Finally, by taking 
the direct sum of Eqs. (7) and (8), a function can be found 
that is only range dependent, and invariant of the orientation 
of the transmitter or receiver: 


(8) 


2 
Vr (m2n1) 


Hi 


2 5 —3cos(2bm1)f 1 17 
a fe ee 
+ VF (n2-2n2) n 


1 2 (9) 
Fy = Fi + Fro =O 
r 


where b=5c? is a constant determined through calibra- 
tions. In practice, F,. is obtained by summing the signal 
power measured at both receivers due to both transmitters. 
The model for inverting range involves an optimization for 
ground conductivity as described in Ref. [3] to obtain a 
model of the range dependence of the geometric distance 
function (1/r°+1/r;), which can be numerically inverted. 
Practical implementation will be described in the following 
sections. 

To obtain a function that is only dependent on orientation, 
the ratio of F,,,. to F,,, can be taken: 


Fr _ 5 —3c08(2hm1) 
Fi 5+ 3c08(2b m1)" 


10 
ae (10) 


where the range dependence is automatically removed 
through cancellation. Solving Eq. (10) for the orientation of 


the transmitter as a function of F4,,,, it is possible to obtain: 


= 1 Tx 
Omi = 3008 


5 -5X Fam 
3X Fam +3 


ay 


1 


> 


where ,,. can be easily obtained if needed by noting that 
9 m2 +7/2. In practice, Fy, is obtained by first summing 
the power of both tones in one receiver, and dividing this 
value by the sum of power of both tones in the second 
receiver. This implementation will be described in the fol- 
lowing sections. 

FIG. 2 illustrates a two-axis transmitter (205) and receiver 
(210) system during range and orientation experiments on a 
building rooftop. Inserts (215, 220) show the related system 
block diagrams. 

To study the accuracy of the decoupled range and orien- 
tation solutions, both transmitter and receiver were devel- 
oped to have two orthogonal axes, and the range and 
orientation of the mobile transmitter was varied relative to 
the fixed receiver system. The experiment was conducted on 
the rooftop of a long building at the Jet Propulsion Labo- 
ratory campus. The setup is shown in FIG. 2. A tape measure 
and orientation table were used to provide range and orien- 
tation ground truth for error analysis. In this exemplary 
system, each axis of the transmitter and receiver consisted of 
about 50-turns of 28 AWG magnet wire wound in a closely 
spaced configuration around a square outline with a side of 
about 20 cm. Both transmit and receive coils were co- 
located at a height of about 1.5 m. The inserts of FIG. 2 (215, 
220) shows the transmitter and receiver block diagrams. 
Two direct digital synthesized (DDS) sources are used to 
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generate tones at 198 kHz and 232 kHz each. The output 
tone of each DDS is amplified using a class-AB amplifier 
design consisting of a LM386 audio amplifier chip-set from 
Texas Instruments (TI). The output of each amplifier drives 
a resonant capacitor network and coil designed to resonate 
at the center frequency of the tones. The emitted magnetic 
field of transmitter coil-1 at about 198 kHz, and of coil-2 at 
about 232 kHz, induces a current in both coils of the 
receiver. To permit detection of both tones in each coil and 
to maintain simplicity in the RF design, a digitally controlled 
switch circuitry was used to switch between coils of the 
receiver. The output of the switch was amplified using a 
two-stage Norton amplifier (Z10043A and Z10047A) from 
Clifton Laboratories. The amplified signals were digitized 
using a +0.1 V, 16-bit 2MS/s analog-to-digital convertor 
(ADC) within the USB-6361 X-Series Data Acquisition 
system from National Instruments. The USB-6361 X-Series 
Data Acquisition was also used to set the digital state of the 
RF switch. Data was streamed to a PC and processed in 
Matlab by Mathworks. Unknown gain/losses in each chan- 
nel of the transmitter and receiver were measured through 
coupling measurements at very short ranges of about 2 m, to 
reduce effects of secondary fields that would distort the 
calibrations. 

The digital time-domain signals of about 0.1 s in length 
captured by the ADC were sent to the PC and converted 
using a Fast-Fourier Transform (FFT) algorithm into the 
frequency domain. The peak power found close to 198 kHz 
and 232 kHz were captured and calibrated using gain 
measurements described previously. By selecting the switch 
digitally, the tones were captured in alternating fashion 
between both receive coils. The measurements during one 
round contained power of each tone at both receiver coils. 
The captured data from the two receiver coils were distin- 
guished by digital switch selection. Since power varies 
linearly with the voltage squared, all 2x2 or 4 power 
measurements were summed within each round. The total 
summed power has the form of F,. in Eq. (9), but requires a 
calibration to remove the unknown gain, b. This is done 
through a measurement at 18.3 m (60 ft). To solve for the 
ground conductivity and a within ra measurement of F,, was 
repeated at 10 m and used to minimize for a, from which a 
ground conductivity of 0.45 S/m was found for the rooftop. 
Next, a model of the geometric range function squared 
((i/P+1/r7 7?) was obtained numerically as a function of 
radial distance, p. The measured total power after calibrating 
F, was then inverted based on this model to obtain the radial 
range, 0, between the transmitter and receiver. 

To obtain the orientation of the transmitter, the total 
summed power of both tones at one receiver was divided by 
the total summed power of both tones at the second receiver. 
This ratio had the exact form of F,,, in Eq. (10) and did not 
require any further calibration. Using Eq. (11), the ratio 
given by F,,, was converted into the orientation of the 
transmitter. 

The above system and methodology are intended as an 
example. Other systems with different configurations may 
also be employed. 

FIG. 3 illustrates measurements of the orientation-invari- 
ant range function compared to theory. The resulting 
inverted range is also compared to ground truth. As visible 
in FIG. 3, the theoretical prediction and experimental mea- 
surements are very close. 

To study the orientation-invariant ranging accuracy, the 
radial range, p, of the transmitter was varied from 12.2 m (40 
ft) to 42.7 m (140 ft) at steps of 10 ft (3.048 m), while 
attempting to maintain an orientation of 9,,,=90°. The total 
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power of the orientation-invariant range function, F,, is 
shown in FIG. 3 (305) to match well with the geometric 
range function squared ([1/r°+1/r;7|*) from complex image 
theory. The solution is numerically inverted to obtain the 
range, which is compared to the expected ground truth range 
in FIG. 3 (310). 

Range error and orientation inversion for measurements 
with variation in range but with fixed orientation are illus- 
trated in FIG. 4. The error in range estimation shown in FIG. 
4 (405) was below 0.2 m for distances up to about over 40 
m. The orientation solutions for the range measurements 
described here varied by about 7° from the expected value 
of 90° , as visible in FIG. 4 (410). 

FIG. 5 illustrates measurements of the range-invariant 
orientation function compared to theory. The resulting 
inverted orientation is also compared to ground truth. To 
study range-invariant orientation accuracy, the orientation, 
9,1, of the transmitter was varied from 0° to 90° at steps of 
10°, at a range of 18.3 m (60 ft). The measured range- 
invariant orientation ratio, Fy,,, is shown in FIG. 5 (505) to 
match well with theory. The solution for the orientation, $,,,), 
is found using Eq. (11), which is compared to the expected 
ground truth orientation in FIG. 5 (510). The orientation 
error, shown in FIG. 6 (605), was below 1° for all orienta- 
tions. The range solutions for the orientation measurements 
described here varied by up to 0.4 m from the expected value 
of 18.3 m (610). 

From the experiment and calculations described above, a 
long-range magnetoquasistatic positioning concept was 
developed that permits the decoupling of range and orien- 
tation in the coupling theory. Therefore, the technique 
described above using magnetoquasistatic fields has been 
shown to enable accurate ranging in non-line-of-sight envi- 
ronments. The system, in the embodiments above, was 
optimized for measurements close to the 2D xy-plane. 
Measurements showed peak orientation-invariant ranging 
errors of 0.2-0.4 m for a range of up to over 40 m. The peak 
range-invariant orientation errors found were 1-7°. 

In the following, the technique described above using 
magnetoquasistatic fields is expanded to include decoupling 
of azimuthal direction of a mobile device. The results are 
azimuthal direction sensing for long-range magnetoquasis- 
tatic systems with a peak error of about 6°, independent and 
invariant of the range or orientation of the mobile device. In 
the following, two-dimensional (2D) positioning is demon- 
strated by combining the decoupled range and azimuth- 
direction solutions. 

As discussed above, the non-line-of-sight (NLoS) posi- 
tioning problem is important in both civilian and military 
applications. Existing techniques that utilize high-frequency 
propagating radio waves suffer in NLoS use (for example, 
see Refs. [4, 6-7]). 

Magnetoquasistatic (MQS) fields are weakly perturbed by 
lossy-dielectrics, and thus beneficial when used in NLoS 
positioning. Direction finding systems require to be progres- 
sively moved closer to the transmitter, and therefore cannot 
provide long range positioning, see Ref. [10]. Short-range 
MQS systems use a frequency of a few kilohertz to avoid 
ground effects, see Refs. [8, 9]. However, the use of very low 
frequency and increased ground effects with range, limits 
their maximum operation range to a few meters. To enable 
long-range NLoS sensing, a MQS technique was developed 
to account for ground effects and by using a higher-fre- 
quency (300-400 kHz) to obtain a range of up to about 50 m, 
see Ref. [3]. The long-range MQS technique enabled rang- 
ing with errors of 0.12-0.7 m for ranges up to 40 m, see Ref 
[3], and was not sensitive to the LoS blockage, see Ref. [1]. 
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However, due to the non-linear position-orientation cou- 
pling and ground effects, the technique can have some 
limitations when extended to 2D or 3D, resulting in errors 
which are dependent on positions, orientations, and the 
non-linear solution convergence, see Ref. [2]. 


As described above in the present disclosure, and as seen 
from Ref [11], a long-range MQS technique has been 
described to decouple the range and orientation of the device 
to be located. The technique can account for ground effects, 
obtaining decoupled solutions with low error over long 
ranges. This section extends the technique to provide 
decoupled azimuth-direction sensing, which enables 2D 
decoupled positioning of the MQS mobile device using a 
single receiver system with orthogonal co-located coils. 


FIG. 7 illustrates transmit-receive coil coupling. The 
secondary fields due to the eddy-currents are modeled using 
complex image theory, see Ref. [3]. 


As described above and in Ref. [11], the voltage measured 
at a receiving coil due to a quasistatic magnetic source above 


ground is due to the source (B.) and image (B,) fields as 
shown in FIG. 7, and is given by Eq. (2). In FIG. 7, the 
system comprises a mobile transmit coil (705) and a fixed 
receive coil (710). 


As detailed above, the total voltage expression of Eq. (2) 
can be simplified to give Eq. (4), where z is the receiver 
height above the xy-plane formed at the origin of the 
transmitter, p is the cylindrical radial range, and 6 is the skin 
depth of penetration in the ground. The key assumption 
required to derive (2) is that at short ranges only the source 
field is significant, but p/r~1, whereas at long ranges, both 
source and image fields are important, but p/r,—1. 


To enable 2D positioning using a single receiver and 
transmitter system, a two-axis orthogonal transmit and 
receive system are employed, as shown in FIG. 8. FIG. 8 
illustrates a top (501) and side (802) view illustration the 
two-axis transmitter and receiver system. The transmitter is 
permitted to have an arbitrary azimuth-direction (@,,,),) 
relative to a reference azimuth-angle. The two-axis trans- 
mitter has horizontal magnetic moments denoted by m, 
(805) and m,, (810) whereas the two-axis receiver has 
horizontal magnetic moments denoted by f, (815) and A, 
(820). The orientation of the transmitter, o,,, (825), is noted 
to be relative to the radial vector direction, p (830). The 
azimuth-direction of the transmitter is defined relative to the 
azimuth angle of the receiver f,, and the angle is given by 
,,, as shown in FIG. 8. 


It can be noted that Eq. (4) exhibits a symmetry between 
the transmitter and receiver moments, th and n. To derive the 
decoupled transmitter orientation, as described above and in 
Ref. [11], a transmitter-centric approach was used, where the 
transmitter orientation, o,,,, was found from the ratio of 
Fg F ma! Ft? 


5-5x Fina. ey VF ¢n2-onk) (12) 
4 pod Fy | Fm k 
ml = 5 ; = ; 
2 ay fmt gf Fm Z VFnt-onk) 
Fn . 


Here, the squared voltage from transmitter 2 measured at 
both receivers (k=1,2) is summed to give F,,..=V7¢m2-2n y+ 
Vigo and divided by similar contributions from trans- 
mitter 1, Fea —¥ neaictnon tata 
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To approach a solution for the azimuth-direction of the 
transmitter, ,,, a receiver-centric solution is derived 
instead. Simplifying Eq. (4) for p-th,=cos 9,1. Pf, =Cos $,,, 
and rh, -f, =cos(9,,,;—9,,; =COs ,,, COS o,,,°Sin ,,,;, the volt- 
age at receiver 1 is found, due to contributions from trans- 
mitter 1, and by using similar steps, from transmitter 2: 


4 (13) 
rr 


Ofoa : : 
Vrantem2) = ~ IG Reosdni sind — COSbym1 SiNGy 1 )] X zt 


where $,,9=9 »+7/2 was used to derive V 74,1 — 2) to keep 
1 in both equations. Taking the sum of squares of Eq. (13), 
the receiver-centric version of the solutions at receiver 1 can 
be found: 


1 2 
Fu = pyogenes 
Al 


where c=-jji,c/4z is a constant. In contrast to F,,,, in Ref 
[11] and in the previous embodiment described above in the 
present disclosure, F,,, is only a function of receiver 1’s 
orientation, and is invariant of the transmitter orientations. 
Following similar steps, the contributions at receiver 2 can 
be found: 


ra={ 


The solutions in Eqs. (14) and (15) are the receiver-centric 
coupled equations which are invariant of the transmitter 
orientations, but dependent on range. To further remove 
range dependence, we take the ratio of Eq. (15) to Eq. (14) 
to find the receiver-centric counterpart to Fy,, described 
above for a previous embodiment, as well as in Ref. [11]: 


Viintemt) (14) 


2 5 + 3cos(2¢n1) EE 
2 rs 


2 
+Viintecm2y 


Vitn2emt) (15) 


pate i 
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2 5 -—3co0s(2¢p1) 
Fy 5 +3co0s(26,1) 


(16) 


Through a direct inversion of Eq. (16), it is possible to 
find the orientation of receiver 1, »,,,, which is the azimuth- 
direction of the transmitter if the receiver is in a fixed 
non-rotating configuration: 


5-5 Fy, 
3X Fon +3 


(17) 


> 


_ 1 -l 
dna = 008 


To solve for the azimuth-direction of the transmitter using 
Eq. (17), it is necessary to compute Fy,=F,,2/F,,,. In practice, 
F,,, and F,,, is obtained by summing the calibrated power due 
to both transmitters at receiver 1 and 2, respectively. 

Once the decoupled azimuth-direction of the transmitter is 
found, its 2D location is known in cylindrical coordinates 
from the azimuth-direction combined with the decoupled 
range solutions, p, as described above in the present disclo- 
sure, and in Ref. [11]. The inverted location of the trans- 
mitter, relative to the receiver location, can also be readily 
converted into Cartesian coordinates using the coordinate 
conversions x=pcos 9,,, and y=psin 9,,,. 
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To study the accuracy of the decoupled azimuth-direction 
sensing, a two-axis orthogonal transmit and receive system 
was developed as shown in FIG. 9, with rotation-tables at 
their bases to allow for rotations. At the transmitter (905), 
two direct digital synthesized (DDS) sources were used with 
amplifiers, and fed to resonant coils for increased signal-to- 
noise. The field transmitted by each axis of the transmitter 
was a tone at 198 kHz and 232 kHz, respectively. Both tones 
are detected at each receiver (910) coil. A switch is used to 
selectively sense the signal at the receiver coils, which is 
then amplified and digitized. The setup is described more in 
details in the present disclosure, above. 

To simplify the experiments for azimuth-direction sensing 
of the transmitter, the transmitter’s orientation, ,,,, its radial 
range, p, and its location a, are fixed, and the orientation of 
the receiver, 9, 1s varied. As described previously and 
illustrated in FIG. 8, the orientation o,,, describes the angle 
between p and fi,, which also describes the azimuth-direc- 
tion of the transmitter relative to the azimuth-reference (9 
when aligned to f, direction) if the receiver’s position and 
orientation is fixed. To characterize the technique and inver- 
sion results, some measurements were taken and compared 
to the theoretical predictions. 

First, the azimuth-direction sensing accuracy for the 
transmitter at a fix range, p, and orientation, ,,,, were 
studied. The transmitter was placed at a radial range of 10 m 
from the receiver system. The locations of both the trans- 
mitter and receiver were fixed, and the transmitter’s orien- 
tation, >,,;, Was set to 0°. To simulate the azimuth-direction, 
the receiver’s orientation, ,,,, was rotated from 0° to 90° at 
steps of 10° . The measured F,,, at each measurement step is 
obtained by summing the calibrated measured power from 
both transmitters at receiver 2, and dividing this by the 
summed powers of both transmitters measured at receiver 1, 
as described by Eqs. (14)-(16). The measured result for F,, 
is compared to the expected theoretical value in FIG. 10 
(1005). FIG. 10 illustrates measurements of the azimuth- 
direction function compared to theory for variation in azi- 
muth-direction of the transmitter (1005). The inverted azi- 
muth-direction is also compared to the ground truth (1010), 
with the related error (1015). 

Inverting the measured F,,, for azimuth-direction, 9,,,, 
resulted in close agreement with the expected direction 
referenced through the receiver rotation table (1010). The 
difference between the expected azimuth-direction from the 
rotation table and the inverted solution is the apparent 
azimuth-direction error. This is shown in (1015) as a func- 
tion of the azimuth-direction, where the error is typically 
less than 2°, except when close to the 0° azimuth-direction. 

For this first experiment, the decoupled range inversion 
was also examined using the technique described above in 
the present disclosure as well as in Ref [11] to confirm that 
it remains invariant of the azimuth-direction. This azimuth- 
direction invariant ranging is important to permit a fully 
decoupled 2D solution. From the same data collected, the 
range function, F,, can be calculated, as described in the 
present disclosure as well as in Ref [11]. A single calibration 
measurement was taken at 60 ft. (18.29 m) to remove 
unknown gains in the system, as described in the present 
disclosure as well as in Ref. [11]. The result for the measured 
F,. compared to theory is shown in FIG. 11 (1105) as a 
function of the variation in azimuth-direction. The inverted 
range using the technique described above in the present 
disclosure as well as in Ref. [11], is found and shown in 
(1110), while (1115) gives the error in range. The range error 
is found to be less than about 14 cm. 
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For the second experiment, the accuracy of the technique 
at a few different ranges was studied by locating the fixed- 
orientation (fixed 9,,,) transmitter at distances, p, of 10 m, 
20 m, and 30 m. The azimuth-direction is varied from 0° to 
90° at steps of 10° for a transmitter range of 10 m, 20 m, and 
30 m. FIG. 12 (1205) shows the resulting azimuth-direction 
error found using the same approach described above. The 
invariance of the ranging solution was also examined at the 
various ranges. The result for this range error is shown in 
(1210). The peak azimuth-direction error found was 
approximately 6°, whereas the peak range error found was 
about 17 cm. 

For the third experiment, the inversions due to various 
transmitter orientations, ,,,, was studied instead. In this 
experiment, the transmitter range, p, was fixed at 20 m, but 
the experiment was repeated for transmit orientations, 9,,,, 
of 0°, 30°, 45°, 60°, and 90°. At each transmit orientation, 
1, the azimuth-direction, 9,,,, was varied from 0° to 90° at 
steps of 10°. FIG. 13 (1305) shows the result for azimuth- 
direction error as a function of azimuth-direction for the 
various transmitter orientations. The azimuth-direction error 
is shown to be below about 6° in all cases measured. FIG. 
13 (1310) shows the result for range error as a function of 
azimuth-direction for the various transmitter orientations at 
20 m. The peak range error was about 0.3 m. 

To provide a summary of the experimental results 
described above, with application to 2D location sensing, 
FIG. 14 shows the experimental results summarized in 
Cartesian coordinates. The receiver is located at the origin 
(1405) of the coordinate system. Circles (1410) denote the 
ground truth of the measurements. For measurements con- 
ducted at a range of 10 m and 30 m, the transmitter 
orientation, >,,,,, was fixed at 0°. Measurements conducted at 
a range of 20 m included transmitter orientations, },,,,, of 0°, 
30°, 45°, 60°, and 90°, as previously described. 

As described above, in some embodiments the long-range 
MQS positioning technique is extended to provide 
decoupled azimuth-direction sensing that is invariant of the 
range or orientation of the transmitter, with a peak error of 
6°. The range inversions were invariant of the azimuth- 
direction and orientation of the transmitter, with a peak error 
of 0.3 m. The technique enables decoupled 2D MQS posi- 
tioning. 

In the following section, measurements of the magneto- 
quasistatic coupling between a co-polarized vertical emit- 
ting and receiving coil, both located above earth, are com- 
pared to calculations using complex image theory and the 
exact integral solutions along a direction perpendicular to 
the surface normal axes of the loops. Inverting the theoreti- 
cal expressions of the exact integral solution to estimate 
distance from measured fields resulted in an root mean 
square (RMS) error of 0.089 m and 0.079 m, for distances 
between 1.3-29.4 m and 1.3-13.5 m, which corresponds to a 
21.9% and 28.4% reduction in RMS distance error, respec- 
tively, over complex image theory. 

In the following, the accuracy of the magnetoquasistatic 
technique described above is studied using the exact integral 
solution to describe the exact fields of a finite loop above the 
earth. The magnetic field of the emitting current loop in 
proximity to a conducting earth can be decomposed into the 
primary fields of the loop and the secondary fields of the 
induced eddy-currents in the earth, see Ref. [3]. The total 
field outside the ground can be approximated by: 


A(x, y, D=Hyx, y, z-h)+H (x, y, -z-h-8(1-/)), (18) 


where c denotes the total fields due to use of complex 
image theory, 8=V2],0@ is the skin depth, l, is the per- 
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meability of free-space, a is the bulk conductivity of the 
earth near the source, and co is the angular frequency of the 
field emitted by the source. 

FIG. 15 illustrates a distance experiment setup and mea- 
surement block diagram (1505). Comparison (1510) of 
complex image theory and the exact integral solutions to 
measurements. The difference in power is shown in (1515). 

Due to the use of a vertical emitting loop depicted in FIG. 
15 (1505), the primary source fields and the secondary 
eddy-current fields are co-polarized in Eq. (18), and thus 
summed. H, is the field of the loop, which contains the 
lowest-order dipole field and higher-order loop fields (qua- 
drupole, et cetera), and can be written in spherical coordi- 
nates using spherical harmonics: 


2 (19) 


—~ FS . , i 
Hy = Fy) 2FUM)Pam st OPP — Film + 1)Pins 1 OD, 
m=0 


where | is the driven current in the emitting loop, | is the 


radius of the emitting loop, r=¥x?+y7+z? is the absolute 
distance from the loop to the observation point, m is the 
number of loop terms with m=0 representing dipole fields, 
P is the Legendre polynomial, ®=cos >, and @ is the angle 
that the surface normal of the loop makes from the x-axis. 
The function FC) is given by 


corens pe (20) 


Fig) = er] 


r 


Eq. (19) is valid for r>l, i-e., outside the loop, which is the 
region of interest. For r>>l, only the m=0 term is important, 
and the fields are that of a magnetic dipole, as seen Ref [3]. 
The theoretical description of the voltage at the terminals of 
the receiver shown in FIG. 15 (1505) is found using Fara- 
day’s law and is given by: 


VI =-jopiolt-H Ja, (21) 


where fi is the unit vector normal to the receive loop and 
a is the surface area of the receive loop. Due to the use of 
vertical emitting and receiving loops which are co-polarized 
with surface normals perpendicular to the f direction, the 
radial component of the field in Eq. (19) does not contribute 


to i-H, or the theoretical description of the voltage, V.”. 

The exact integral solutions for the fields of a finite loop 
above earth accounts for all higher-order loop terms (all m) 
and the exact induced eddy-currents. Because complex 
image theory only approximately models the fields of the 
exact induced eddy-currents, it is expected that there will be 
some improvements in distance measurement using the 
exact integral solutions. The magnetic fields using the exact 
integral solution of a finite circular current carrying loop 
located above earth is given in Cartesian coordinate by (see 
Ref. [13]): 


(22) 


a 


il , Bo ‘ 
H, = = f f F.|[ Tax) (py) + — 5 + T(Byt dod B, 
n Ay Ay 


where the subscript e denotes the exact fields of the 
emitting loop above earth, a and £ are integration variables 


which operate on each integrand from 0 to +0, 4,=Va7+p, 
T(C)=tan(C) and F, is given by: 
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where I, is the modified Bessel function of the first kind, 
2», Vo2+B74jx?, and K=V2/8. The theoretical description of 
the voltage at the terminals of the receiving loop is found 
using Faraday’s law: 


VF =joug\i-H,Ja. (24) 


Due to the use of vertical loops with surface normals 
perpendicular to the x direction, the x and Z component of 


the exact field in Eq. (22) does not contribute to 4-H, or the 
theoretical description of the exact voltage, V,”. 

Only the one-dimensional (1D) case shown in FIG. 15 is 
studied in the following, where distance is varied along a 
direction perpendicular to the surface normal axes of the 
loops, i.e., x coordinate is varied. The experiment was 
conducted above earth on an open field and has been 
previously described in Ref [3]. The result for the measured 
voltage is plotted in FIG. 15 (1510) as measured power in 
dB, 10 log,,!W“17/R. The complex image theory description 
in Eq. (21) for dipole fields (m=0), and the exact integral 
description in Eq. (24) are also plotted in (1510), where a 
fixed offset was used to calibrate out the measurement gain 
and to make both theories agree with the measurements at a 
distance of 20.2 m. FIG. 15 (1515) shows the difference 
between the complex image and exact descriptions from the 
measurements. 

FIG. 16 illustrates a comparison of: (1605) estimated 
distance error from inverting the theoretical expressions for 
a magnetic dipole (m=0) and a finite-sized loop (m=0,1,2) 
both using complex image theory, and for the exact integral 
solutions; (1610) RMS distance error reduction using exact 
integral solutions. 

The closed-form field equations in Eq. (19) describing 
complex image theory allows a direct numerical inversion 
for the distance. Here, the emitter’s x coordinate can be 
found using a non-linear least square optimization algorithm 
to minimize [V,’-V™]?, where V” is the measured voltage 
at the terminals of the receiving loop, see Ref. [3]. FIG. 16 
(1605) shows the result of this optimization, where the 
distance error is estimated by finding the difference between 
the x coordinate solutions and the actual x coordinate. The 
black line (m=0) denotes the distance error due to use of 
complex image theory with dipole fields, whereas the blue 
line (m=0,1,2) includes the next two higher-order fields of 
the finite sized loop, see Ref. [12]. 

Use of the exact integral solution, which does not have a 
closed-form, precludes the possibility of inversion as done 
with complex image theory. Instead, by using the exact 
integral solutions in Eq. (24) of a finite-sized vertical loop, 
a 1D look-up table was created with the exact theoretical 
induced voltage values and interpolated to minimize [V,7— 
V™/? and find an estimated x coordinate. FIG. 16 (1605) also 
shows the result of this interpolation and minimization, 
where the distance error is again estimated by finding the 
difference between the x coordinate solutions and the actual 
x coordinate. The line (1615) denotes the distance error due 
to use of the exact integral solutions. FIG. 16 (1610) 
provides an analysis of the RMS distance error using com- 
plex image theory and the exact integral solutions. The exact 
integral solution gives an RMS error of 0.089 m and 0.079 
m, for distances between 1.3-29.4 m and 1.3-13.5 m, respec- 
tively. This is compared to the result of 0.114 m and 0.110 
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m, for the same distances from use of complex image theory 
and spherical harmonics to model a finite-sized loop with 
m=0,..., 4. 

As discussed above, the magnetoquasistatic field equa- 
tions can be decoupled into position and orientation in the 
two-dimensional (2D) space of the fixed sensor or mobile 
device plane, even over long-ranges (see Ref. [11]). This 
approach is possible due to a few mathematical approxima- 
tions made possible by the physical nature of magnetoqua- 
sistatic coupling near a conducting ground (as described 
above and in Ref. [11]). The device to be tracked can be the 
transmitter or the receiver, with the fixed sensor being the 
receiver or the transmitter. In addition to the 2D position and 
orientation of the device, the theory can be further extended 
to enable 2D azimuth direction sensing of the device (as 
discussed above and in Ref [14]). A remaining problem of 
importance is the case of z-location tracking in the cartesian 
3D coordinate frame. Under this translation, it is important 
to know if the magnetoquasistatic field equations can be 
decoupled meaningfully for 3D position sensing. 

In the following section, the 3D problem in magnetoqua- 
sistatics is discussed and it is shown that the field equations 
can be decoupled to enable linear 3D position sensing of the 
device by extending the framework presented above. The 
closed-form solutions derived here permit full 3D position 
sensing of the mobile device using a single fixed device, and 
do not require trilateration or triangulation. 

FIG. 17 illustrates the 3D problem with 3-axis transmit- 
ting (1705) and receiving (1710) coils and their respective 
magnetic moments. The complex image is also shown below 
the ground (1715). 

FIG. 18 illustrates a mathematical description of the 3D 
position sensing problem. The receiver (1810) is assumed as 
the device to be tracked, without loss of generality. The 
device to be tracked is permitted a change in z-location. A 
transmitter (1805) is also illustrated, with a ground image 
(1815) and a top view (1820). In contrast with FIG. 1, the 
magnetic moments are illustrated with three components 
instead of two. 

In the previous embodiments, two-axis coils were used. 
However, in other embodiments, three-axis coils can be 
used, permitting 3D instead of 2D positioning. 

For 3D position sensing in MQS over long-ranges, ground 
effects due to induced eddy-currents in the ground cannot be 
neglected. Following the framework described above and in 
Ref. [11], it is required that the transmitting and receiving 
device have 3-axis coils. The system is described in FIG. 17, 
where each transmitter coil and its corresponding complex 
image is described using magnetic moments to simplify the 
mathematical discussion. This approximation 1s sufficiently 
accurate for any appreciable range (greater than a few 
diameters from the coils). 

To understand the effect of 3-axis magnetic moments, it is 
necessary to study the coupling at the receiver due to an 
arbitrarily oriented transmitter magnetic moment. Due to the 
boundary condition made by the air-ground interface, it is 
helpful to decompose the arbitrary transmitter moment into 
its vertical and horizontal components, m=m,p(>)+m,z, 
where the magnitude of the moment in the horizontal plane 
is given by m,, while the magnitude of the moment in the 
vertical plane is given by m,. Following this convention, it 
is also helpful to decompose the receiver’s moment into its 
vertical and horizontal components, A=n,p(p)+n.2. The volt- 
age measured at the terminals of the receiver coil with the 
surface normal, f, is given by Eq. (25), where c=-jap,a/4a. 
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The terms in Eq. (25) for the horizontal and vertical 
moments do not act on perpendicular f unit vectors, there- 
fore Eq. (25) can be significantly simplified. Additionally, 
since the vertical magnetic moments of m,z interact in a 
destructive manner with its complex image, —m,1m, it is 
useful to gather these terms separately from all other terms. 
Finally, to simplify the analysis, two separable cases are 
studied below: the total voltage at the horizontal moment 
sensor, n,O, and at the vertical moment sensor, n,~. With 
these voltages, the total voltage sensed at the horizontal II is 
given by: 


Vip) = (26) 


(e/r? 


Pp 


(p/ ri) 
ia a 


OMo®!, _, , rap tos. 1 oll 1 
lao, ‘Pip | |- Gi, -tip| + =| + 


pel?) ipl?) 
rn a 


3(mi,- 2)(Tp | 


The third term in Eq. (26) is negligible for medium to long 
ranges or distances, such that (r/r,°—>1, if the height is 
approximately less than the radial distances, such that z<p. 
For longer ranges or distances, the third term can also be 
neglected for p<z. The main reason for this is that since a 
3-axis transmitter moment will be used, the first term, due to 
the horizontal moment, will have a much stronger contri- 
bution than the third term due to the cancelation in the third 
term. This is always true since the two horizontal moments 
in the transmitter are orthogonal to each other and will 
always give a strong contribution to the first term. Under 
these criteria, the total voltage sensed at the horizontal II can 
be simplified to give: 


V(ftp) = (27) 


Ope, F(pir)? (piri 
42" lao, Pip pee + ia. = 


i 


ee hfe Sa 
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The goal as discussed in Ref. [11] is to approximately 
decouple the range terms from the orientation, elevation, and 
azimuth terms in Eq. (27). To this order, Eq. (27) can be 
rewritten by a simple factorization as: 


P 5) ot bao a Palas 1 (28) 
Vifip) a cB Mp ‘Pp “PCy = (Mp-Np)] Pp a =| 
where c=jop,0.4/r, and the factor Cy is given by: 
_(elrh [rr tein IF (29) 
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where it follows that Cy~(p/ry? is a good approximation 
for short to long ranges. In the following, the remaining 
vertical equations are discussed. 

The total voltage sensed at the vertical fi is given by: 


o, (30) 


2 
6a, pre] 2 a 
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Following similar considerations (and as can be seen in 
Ref [14]), it can be noticed that under any appreciably 
long-ranges, the third term becomes negligible. The second 
term can also be neglected for appreciable ranges since it 
will be much weaker than the second term in Eq. (27). With 
these approximations, the total voltage sensed at the vertical 
fi can be simplified to give: 


(p2/) | (2/77) G1) 


rn r 


Vii.) = -j-E [Gri Pi, -2)] 


It can be noticed that Eq. (31) can be rewritten by simple 
factorization to obtain: 


ee eas 32 
V(Az) = c[3Mip- pP), -2)Gi| 5 = =| me 
where the factor C, is given by: 
cy = LPIA + 0%, 1DIP (3) 
as 1/P+1/7A 
_ (pz/P {Pe +0 _ pz 
<"T7R+0 Pe 
2x(pz/n2/P_ pz 
y oxi = roe row, 


where it follows that Cg=pz/r? is a good approximation for 
short to long ranges. The above concludes the derivations for 
the MQS field theory. The solutions derived in Eqs. (28) and 
(32) are the foundations of the following section, where the 
positioning framework is presented. 

The orthogonality relations of the transmitter moments 
are always preserved under any arbitrary rotations, due to 
the use of orthogonal coils on the transmitter, orthogonal 
magnetic moments m,, mh, and m,, and field-magnitude 
only sensing of in the receive coils. As discussed above, it 
can be assumed that m,, fil,, and in, are all mathematically 
similar, or in other words that these quantities have the same 
gains (or calibrated gains) as well as close center frequen- 
cies. With this condition, the preservation of orthogonality 
implies that under any arbitrary rotation there will be two 
horizontal magnetic moments, which are orthogonal to each 
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other and parallel to the horizontal plane, and one vertical 
magnetic moment parallel to the vertical plane. It is then 
possible to obtain the solutions to the elevation-angle direc- 
tion of the transmitter, and then the azimuth-angle direction 
of the transmitter. In a first step, the transmitter and receiver 
centric F equations are defined for the 3D case, which is 
related to the 2D case, as described above and in Refs. 
[11,14], where x=1,2,3, referring to each transmitter 
moment: 


Prax=Vonx—-n pt Vmx—=n2 + Vowssteay? (24) 


= 2 2 2 
Fax=V oax—mi) +Vax—miy + VM Gucemx) + 


For 3D positioning, in some embodiments it is possible to 
focus on the receiver-centric F,. Under any rotations, as 
described above, there will be two orthogonal horizontal 
transmitter magnetic moments. These moments can be 
denoted as located along m, and m,, which are orthogonal 
to each other, but can assume any arbitrary orientation from 
P. Some orientations for this problem can be defined as 
follows (and also as in FIG. 17): 


My P = COSPmi (35) 
fy Pp = cosdny 
my Ay = cos(bmt — Ont) 


= COSPm1 COSPn! + SIND n1 SIND»! - 
Based on this definition, we find the voltage due to 


transmitter m, at receiver fi, as: 


(36) 


Vint-nt = 


C[3Cgcos bmi COSGn, — COSPm1 COSPn1 — SING» SiN, 1 ] X 


1 1 
anal 


which can be re-written in a more convenient form for 
further manipulations: 


; 3 ime (37) 
Vini-nt = C[ACOSbn1 COSH, 1 — sity sind | a =| 
aa 
where A=3C,-1. Similarly, it is possible to find V,,,5_,1: 
(38) 


4 4 1 ol 
Vino—nt = C[Acosh»2C08hn1 — sinbyesind | = =| 
7 


Since the vertical transmitter moment does not contribute 
to the horizontal receiver moments as indicated by Eq. (28), 
V n3-m1=0. It is then possible to find through the same 
simplification approach shown above and in Refs. [11, 14], 
the receiver-1 -centric equations, F,,,. 


Vont-nty + G9) 
2 
rt 2 one. 
Fr = 4 Vino-ny + p= OBL Ag ; 
Ver3-—nl) 
where B, is given by: 
(A? +1) -(A? = 1)cos26y1 (40) 
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Following the same steps, it is possible to find the 
receiver-2-centric equations, F,,,. 


Vout + (41) 
1 if 
BR 


2 2 
Fro = 4 Ven2-n2) + ¢ = By 


> 


2: 
Vin3—n2) 


where B, is given by: 


ie (A? +1) - (A? - L)cos2¢y1 (42) 


2, 5) > 


where the substitution »,.=o,,,+7/2 was made to permit 
writing F,,, as a function of 9,,,. 

The final equation needed is that for F,,, which is the 
vertical receiver-3-centric equation and can be found using 


Eq. (32). To obtain F,,3, Vj. -n35 Wmo-n3, ad V,,3-73 are 
calculated: 
Lot (43) 
Vint-n3 = C[3Cgcoshynt ] pts 


1 1 
Vind—n3 = c{3Ci€08¢ya| + =| 
Poor 


Vina—n3 = 0. 


Using the relation 9,,.=,,,+7/2, and the identity 
cos”,,+8in7},,;=1, F,,3 is then calculated: 


Vent-n3) + i iss (44) 
Fra =) Vina-n3) + f= |= + =| ; 
2 ied 
Viin3-n3) 
where B, is given by: 
B3=(3Cg). (45) 


Having found all three receiver-centric equations, F,,,, 
F,,,, and F,,,, all critical equations are available to solve for 
the 3D positions in MQS for short to long ranges. 

In a next step the method continues by solving for the 
elevation-direction angle of the transmitter. To accomplish 
this, the ratio between the vertical to horizontal F,, is used. 
In particular, from Eqs. (39), (41), and (44), it can be found 
that: 


Fra BC” (BCP (46) 


FutFy A?+1 BCy-1P 41 


and, as previously found, Cg~pz/r* and Cy~(p/r). Now, 
defining the azimuth-direction angle as cos 8=p/r, it can be 
found find that z/r=sin 0, from which it can be obtained: 


Cgxcos Osin 0 


Cyecos*@, 


such that the ratio F,./(F,,,+F,,.) can be written as: 


Fy (3cos@sin@)” (48) 


FutFo @coe@—1 +1 
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FIG. 19 illustrates a plot of F,,;/(F,,,+F,,2), the elevation- 
direction angle ratio function, showing an ambiguity in the 
solution space beyond about 52° (1905). This is not a major 
concern in most ground based position location systems 
operating over long ranges. 

The next step in the method is to find the azimuth- 
direction angle of the transmitter. This can be found by again 
looking at the receiver-centric F,, equations. In particular, as 
described above and in Ref. [11], it is known that the ratio 
of F,,/F,,. gives a solution for the azimuth-direction angle in 
2D. For the 3D case, the ratio of Eqs. (39) to (41) can be 
calculated to find: 


Fy (49) 
Fy = Fo = 


BL (A? +1)-(A? = 1c0826 1 
BQ (A? +1) +(A? = l)cos2¢,1’ 


which can be directly inverted in a similar fashion to the 
2D case (as described above and in Ref. [II]) to give: 


1 f(A? +)-(74+)xFy (S0) 
oni = =Cos | —s————_——__. |, 
2 (A2 = 1) x Fy + (A? -1) 
where it can be recalled that A=3Cg-1: 


A=3C,-3 cos’0-1, (51) 


and where it can also be noted that, when 0=0, the 
azimuth-direction angle solution in Eq. (50) becomes pre- 
cisely that derived in the 2D case (as described above and in 
Ref [11]). 

The above equations give the solutions for the elevation- 
direction angle and the azimuth-direction angle of the trans- 
mitter. The only remaining item to obtain the full 3D 
location is the range. Following the approach described 
above, and in Ref. [11], it is possible to sum the contribu- 
tions from all combinations of the transmitter moments to 
the receiver moments. This can be obtained by summing the 
contributions from all receiver-centric functions, F,,,+F,,.+ 
F.5: 


5 (52) 
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F,= 4 Fro + —+— 
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2 1 1 
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where the term in the brackets can be rewritten as: 
A? +1 +(3Cg)? = (3Cg - 1)? +1 + (3C9)? (53) 
= (3cos?6 — 1y + (3cos@siné)? + 1 


=5-3sin76, 


which results in the simplified F,.: 


1 if (54) 
st 


F,=c7[5-3sin’@]x| = +5 
ti Fj 


For 0=0, F,, reduces to exactly the solutions derived in the 
2D case, as described above and in Ref [11]. This completes 
the solutions for the 3D position of the transmitter. 

The algorithm to determine the 3D positions is as follows: 
By numerically solving Eq. (48), a solution is obtained for 
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the elevation-direction angle, 8. Once 0 is known, A is also 
known, from which the inverted form in Eq. (50) can be used 
to solve directly for the azimuth-direction angle, and Eq. 
(54) can be used to solve for the range to the transmitter 
through the same numerical inversion technique described 
above and in Ref [11] for the 2D range. 

In the above section, a new decoupled theory in magne- 
toquasistatics was described, showing that the 3D positions 
of a mobile magnetoquasistatic device can be found using 
linear closed-form equations. The linear closed form solu- 
tions can be used to find the 3D location of a mobile device 
over short to long ranges above a half-space such as earth. 

In the following the 3D orientation sensing problem of the 
mobile device is discussed. In the section above the 
decoupled magnetoquasistatic coupling between an arbi- 
trarily oriented mobile magnetic dipole (or an excited elec- 
trically-small coil) and an arbitrary receiving coil was 
derived. Due to the use of a 3-axis receiver architecture, the 
problem was simplified by assuming that the receiving coils 
were either horizontal (vertical moment) or vertical (hori- 
zontal moment). The purpose was to decompose the problem 
in relation to the boundary condition to enable a simplifi- 
cation in the analysis. 

The coupling from the transmitter moment to the hori- 
zontal receiver moment is given by (as seen above and in 
Ref. [15]: 


A er Z Byonss 1 1 (55) 
V(fip) = c[3(Mip ‘Pri “PCy — (Nip “igs + =|. 


where c=jwil,a/4a, and the factor Cy~(e/ry. The cou- 
pling from the transmitter moment to the vertical receiver 
moment is given by: 


(56) 


> Bn ay 1 
V(Ai,) = c[3m, “Piz 2)Co - (Mp -Mp)] Pp + 
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where the factor Co~pz/r’. Eqs. (55) and (56) are general 
to arbitrary orientations of the transmitting magnetic 
moment. To enable closed-form solutions for 3D position 
sensing, it was shown that both the transmit and receive 
system required 3-axis coils. To simplify the analysis, it was 
chosen to constrain the 3-axis transmitter system such that 
its moments were aligned along the x, y, and z axis, 
respectively. This resulted in th,-p=cos ©,,;, and th, “f,=cos 
(911-91); for the magnetic moment of the transmitter 
moment, m,. To remove this restriction, the dot product is 
redefined as follows: 


My +P = COSOm1COSGm1 (57) 


fy -P = Cosby 
7 +My = COSHn1COS(bm1 — Prt) 


= COSA nj [COSGn1 COSGy1 + SING m1 SING! |, 


where the factor cosO,,,, describes the elevation of the 
transmitting magnetic dipole from the horizontal plane. 

Due to Eq. (57), the explicit coupling equations need to be 
updated. To obtain range-invariant orientation sensing, the 
transmitter-centric coupling equations need to be studied, 
see Ref. [11]. The coupling of m, to the each of the 3-axis 
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coils is studied, now with the assumption that m, can obtain 
an arbitrary 3D orientation based on Eq. (57). The voltage 
measured due to m, at receiver n, is given by: 


Vint-nl = (8) 


il 1 
C[3CgcOSAn1 COSPm1 COSGn1 — COSPm1 COS(Pm1 — Hpi) X FE + =| 


Following similar manipulations, the voltage at receivers 
n, and n, can be found: 


Vini-n2 = (59) 


1 1 


C[3CyCOSA;n1 COSPn1 COSHn2 — COSA 1 COS(Pint — Gn2)] X| | + ah 
i ? 


Pe 
=t+s|. 
P rn 


(60) 
Vint-n3 = C[3CeCOSA m1 COSHm1 | X 


The solutions above describe the coupling to each 
receiver moment due to an arbitrarily located and orientated 
transmitting magnetic moment. Therefore, the solutions for 
the remaining magnetic moment, m,, can be found by using 
0,,2 instead of o,,,; and ,,,> instead of @,,,, and similarly, the 
solutions for m,, can be found by using 8,,, instead of 6,,,, 
and ,,3 instead of 9,,; 

The new solutions in Eqs. (58)-(60) are different com- 
pared with those in Ref. [15]. It is expected that the position 
and ranging functions are still invariant under arbitrary 
rotations of the 3-axis transmitter system, due to the ortho- 
normal basis used. It is valuable to show that this is indeed 
the case with Eqs. (58)-(60). The total power due to m, is 
first found, where upon some mathematical manipulations 
and simplifications, it is obtained: 


Ven) + (61) 
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where B,,,, is given by: 
By 1=COS7O,,1 COSPyy +(9CG?—6C 99 +9C 199 7)COS"Pyn + (62) 


Following the same approach to find F,, , for m,, and 
F, 3 for m,, the sum F,=F,,,+F,,.+F,,,3 can be found: 


Font) + : (63) 
flied 
F, = 4 Fema) + ¢ = C° Bm ata : 
Fim) 
where B,,, is now given by: 
Bm = COS?O1 COS bmi + COS” On2C0S 2 bm + COS” On3COS? Pind + (64) 


9c3 —6Cg + 95) X [cos*Gmt + COS? bmg + C08" bn3] = 


1+9C5 -6C,y +95, 


where C,=cos0 and C,g=cos Osin 0. The range solutions 
that are derived from F,, in Eq. (63) (see Ref. [15]) are only 
a function of the elevation-direction angle of the transmitter, 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


65 


24 


@, and not a function of the transmitters orientations, 0,,,1 9,3 
OF ®»1,9,3- 4 similar analysis can be done to show that F,,, 5 3 
in Ref. [15], and hence the azimuth-direction and elevation- 
direction angles remains invariant of the transmitters orien- 
tations, 8,,1.5.3 OF Pmio,3- 

To emphasize the azimuth orientation of the transmitter 
and based on the form of Eqs. (55) and (56), it is possible to 
analyze the ratio of V 21V. ? to find: 


ml-nl ml-n3 


Ves _ (c0SOn1 — 3Cgcosdny + sing, tanbyn1) (65) 
Vink-n3 9C3 ; 


Which is dependent on 9,,,, along with 9,,, and 0. To 
remove the dependence on 9,,,, it is possible to use the 
summation of the power of ratios between the horizontal and 


vertical coupling equations, Vyji—17/Vmniens-tV mi-no! 
2, 
Vini-n3 * 
Vanni, Vata _ 9C4 ~ Cp + Sec” Pmt (66) 
Fam = 2 Veg 9G : 
mi—n3 mi—n3 @ 


which is no longer dependent on 9,,,. Eq. (66) can be 
directly inverted to find the azimuth orientation of the 
transmitter moment m,: 


(67) 


1 1/2 
= cos7! : 
Pmi = COs la =9C3 + 6C, 


To emphasize the elevation orientation of the transmitter 
and based on the form of Eqs. (55) and (56), it is possible to 
investigate the ratio of V ?/F.. to find: 


m1—n3 


Fn = Ve 7 9CZcos* Pn Si” Ont (68) 
oe > “aie 9C} +(1 -3Cg)? * 


where we find dependence on cos 9,,, and sin 8,,,. Since 
it is possible to solve for cos o,,,, using Eq. (67), it is possible 
in-turn to invert for the elevation angle of the transmitter 
moment m,, 9,,,, by inverting Eq. (68): 


12 
Fon | 


The derivations given in Eqs. (67) and (69) are straight- 
forward and elegant, but present a problem. Firstly, Eq. (67) 
is not sensitive to Fy,,,; when C,>0, which is the case on the 
2D x-y plane that is of importance in long-range ground 
based positioning. Secondly, Eq. (69) is not sensitive to F,,,, 
when C,-0 or cos 9,,,,>0, which is also on the 2D x-y plane 
or for $,,,;>7/2. The reason for this difficulty is that the 
1/,,2;__3 is used, which becomes a null field on the 2D x-y 
plane for long-distances due to multiplication with C,=pz/r’, 
caused by the boundary conditions. In addition, it is not 
possible to use the ratio conventions used in Ref [11] 
because the magnetic moment is no longer on a fixed 2D 
plane. To avoid these problems, it is possible to take the ratio 


mi? 


(69) 
Om = sin! 


|= +(1-3Cy)* 


9C2cos bmi 
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between the field seen at receiver coil 1 to receiver coil 2, 
since both have horizontal moments to find: 


Vind—nl (79) 
fal Vint-n2 
_ (©£08Prmt = Pn) = 3CgCOS(Pmt EOSPrt »? 
(3Cgcos(Gnt sin) + sina — Oa)?” 


which is now a function of __,,;, $,,, and Cy. In Eq. (70), 
Cg is not present and does not reduce the sensitivity of F,,,,; 
in the inverted solutions. The un-approximated solution for 
$1 18 then given by: 


ss a D1+D2+D3 
Pmt =—CO8 | Dae DS 


D1 =4(Fon — Vcos(2¢nt) 


(71) 


> 


D2 = 32y Fam (Fat + 1)?sin? (2b) 
D3 = -3(Fami + 1? cos(4Gnt) + Figmt (7 Fam +38) +7 


D4 = -60(F 5m — Lcos(2bn1) 


DS =9(Famt + 1)? cos(4nt) + Femi (S9Fami +46) + 59. 


The closed form solution given in Eq. (71) is always valid 
and sensitive to Fy,,, under any Cg, Cy, and for all ,,,. Once 
9,1 18 known, it is possible to find the elevation angle of the 
transmitter moment by instead studying the ratio of 
Vinten tV mien +V¥ mina to the total power found through 
F.. in Ref. [15]: 


Vint_nt + Vintna + Vintn3 (72) 


Fant = 
Om1L F, 


_ £08" (Bm (3(3C§ + Cp(3Cy — 2))c08" (Fmt) + 1) 
> 9C3 + (1 -3C4 . 


which can be directly inverted to find the elevation angle, 
0 . 


mi* 


V Fam ¥9C3 +9C3 -6Cy +1 ) 


VD6+D7 


On, = COS “| — 


D6 = 9C5cos? (bm ) 


D7 = 9C§c08" (Pm) — 6CgCOs? (Pm) + 1. 


The algorithm is straightforward. By solving Eqs. (71) 
and then Eq. (73), it is possible to find the azimuth orien- 
tation and elevation orientation the transmitter moment m,. 
It is possible to repeat similar calculations to obtain Fy, and 
F4,,2 for transmitter moment m,, from which it is possible to 
solve for 9,,> and 0,2, and then F,,,,, and F,,,3 for transmitter 
moment m,. In turn, it is then possible to solve for ,,, and 
0,,3- With these solutions, the full arbitrary 3D orientation of 
the mobile transmitter system is known. 

The current closed-form solutions for 3D orientation 
sensing, used with the closed-form solutions for 3D position 
sensing derived in Ref. [15], completes the problem of 3D 
position and orientation sensing using magnetoquasistatic 
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fields for short to long ranges using complex image theory. 
The forms derived allow use of complex image theory or 
higher order corrections to complex image theory, as noted 
in Ref [11], and permit closed-form solutions for 3D position 
and orientation sensing of a mobile device using a single 
fixed sensor with 3-axis orthogonal coils. 

In some embodiments of the present disclosure, an advan- 
tage of the methods described herein is the ability to use a 
single sensor instead of multiple sensors, to carry out the 
positioning through quasistatic magnetic fields. Another 
advantage is the decoupling of the equations for the orien- 
tation and the position, therefore having decoupled errors in 
orientation and position. In some embodiments, the methods 
of the present disclosure can be carried out by having a 
receiver or a transmitter with a known location. This can 
mean that the receiver or transmitter have a fixed location, 
or they may be moving in a known way, therefore still 
having a known location. In embodiments with two or three 
axis, implemented, for example, with two or three orthogo- 
nal coils, each coil can sense a signal from the two or three 
coils in the transmitter. In other words, a first receiver coil 
can independently sense two or three independent signals 
from each of the two or three transmitting coils. Several 
methods can be employed to maintain the independence 
between the two or three signals. For example, a different 
frequency may be used for each of the three orthogonal 
coils, or multiplexing may be used. Multiplexing can com- 
prise sending signals for each direction at a different time, 
according to a specified sequence. Other methods to differ- 
entiate between the three signals can also be implemented, 
as understood by the person of ordinary skill in the art. In 
some embodiments, the methods described in the present 
disclosure decouple the horizontal and vertical components 
of the ground fields. 

The methods and systems described in the present dis- 
closure may be implemented in hardware, software, firm- 
ware or any combination thereof. Features described as 
blocks, modules or components may be implemented 
together (e.g., in a logic device such as an integrated logic 
device) or separately (e.g., as separate connected logic 
devices). The software portion of the methods of the present 
disclosure may comprise a computer-readable medium 
which comprises instructions that, when executed, perform, 
at least in part, the described methods. The computer- 
readable medium may comprise, for example, a random 
access memory (RAM) and/or a read-only memory (ROM). 
The instructions may be executed by a processor (e.g., a 
digital signal processor (DSP), an application specific inte- 
grated circuit (ASIC), a field programmable logic array 
(FPGA), a graphic processing unit (GPU) or a general 
purpose GPU). 

A number of embodiments of the disclosure have been 
described. Nevertheless, it will be understood that various 
modifications may be made without departing from the spirit 
and scope of the present disclosure. Accordingly, other 
embodiments are within the scope of the following claims. 

The examples set forth above are provided to those of 
ordinary skill in the art as a complete disclosure and descrip- 
tion of how to make and use the embodiments of the 
disclosure, and are not intended to limit the scope of what 
the inventor/inventors regard as their disclosure. 

Modifications of the above-described modes for carrying 
out the methods and systems herein disclosed that are 
obvious to persons of skill in the art are intended to be within 
the scope of the following claims. All patents and publica- 
tions mentioned in the specification are indicative of the 
levels of skill of those skilled in the art to which the 
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disclosure pertains. All references cited in this disclosure are 
incorporated by reference to the same extent as if each 
reference had been incorporated by reference in its entirety 
individually. 

It is to be understood that the disclosure is not limited to 
particular methods or systems, which can, of course, vary. It 
is also to be understood that the terminology used herein is 
for the purpose of describing particular embodiments only, 
and is not intended to be limiting. As used in this specifi- 
cation and the appended claims, the singular forms “a,” an 
and the include plural referents unless the content clearly 
dictates otherwise. The term “plurality” includes two or 
more referents unless the content clearly dictates otherwise. 
Unless defined otherwise, all technical and scientific terms 
used herein have the same meaning as commonly under- 
stood by one of ordinary skill in the art to which the 
disclosure pertains. 

The references in the present application, shown in the 
reference list below, are incorporated herein by reference in 
their entirety. 
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What is claimed is: 

1. A method comprising: 

providing a receiver, the receiver comprising a first and 

second receiver coils configured to receive quasistatic 
magnetic fields in a first and second mutually orthogo- 
nal directions transmitted by a transmitter, the trans- 
mitter comprising a first and second transmitter coils 
and configured to transmit the quasistatic magnetic 
fields in a third and fourth mutually orthogonal direc- 
tions; 

receiving the quasistatic magnetic fields through the 

receiver; 

measuring, with the receiver, a first voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the first receiver 
coil; 

measuring, with the receiver, a second voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the first 
receiver coil; 

measuring, with the receiver, a third voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a fourth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the second 
receiver coil; and 

based on the first, second, third, and fourth voltages, 

calculating a distance between the receiver and the 
transmitter, wherein the measuring of each of the first, 
second, third, and fourth voltages is independent from 
the remaining voltages. 

2. The method of claim 1, further comprising, based on 
the first, second, third, and fourth voltages, calculating an 
orientation of the receiver or of the transmitter. 

3. The method of claim 1, further comprising, based on 
the first, second, third, and fourth voltages, calculating an 
azimuthal angle of the receiver or of the transmitter. 

4. The method of claim 2, wherein calculating a distance 
between the receiver and the transmitter and calculating an 
orientation of the receiver or of the transmitter are based on 
mutually decoupling the distance and the orientation. 

5. The method of claim 4, wherein providing a receiver 
comprises only one receiver system. 

6. The method of claim 4, wherein calculating a distance 
between the receiver and the transmitter and calculating an 
orientation of the receiver or of the transmitter comprise 
receiving a location of either the receiver or the transmitter. 
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7. The method of claim 6, wherein independently mea- 
suring each of the first, second, third, and fourth voltages is 
by receiving the quasistatic magnetic fields at a different 
frequency for the mutually orthogonal directions. 

8. The method of claim 6, wherein independently mea- 
suring each of the first, second, third, and fourth voltages is 
by receiving the quasistatic magnetic fields at a different 
time for the mutually orthogonal directions. 

9. The method of claim 2, wherein calculating a distance 
between the receiver and the transmitter and calculating an 
orientation of the receiver or of the transmitter are based on 
decomposing a near-ground magnetoquasistatic dipole into 
horizontal and vertical components. 

10. The method of claim 9, wherein decomposing a 
near-ground magnetoquasistatic dipole into horizontal and 
vertical components comprises approximating the near- 
ground magnetoquasistatic dipole by complex image theory. 

11. The method of claim 10, wherein calculating a dis- 
tance between the receiver and the transmitter and calculat- 
ing an orientation of the receiver or of the transmitter are 
further based on decomposing the received quasistatic mag- 
netic fields into horizontal and vertical components, and 
using the horizontal component only, for the received qua- 
sistatic magnetic fields. 

12. A method comprising: 

providing a receiver, the receiver comprising a first and 

second receiver coils configured to receive quasistatic 
magnetic fields in a first and second mutually orthogo- 
nal directions; 

providing a transmitter, the transmitter comprising a first 

and second transmitter coils and configured to transmit 
the quasistatic magnetic fields in a third and fourth 
mutually orthogonal directions; 


transmitting the quasistatic magnetic fields through the 
transmitter; 

receiving the quasistatic magnetic fields through the 
receiver; 


measuring, with the receiver, a first voltage originating 
from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the first receiver 
coil; 

measuring, with the receiver, a second voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a third voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the first 
receiver coil; 

measuring, with the receiver, a fourth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the second 
receiver coil; and 

based on the first, second, third, and fourth voltages, 

calculating a distance between the receiver and the 
transmitter, wherein the measuring of each of the first, 
second, third, and fourth voltages is independent from 
the remaining voltages. 

13. The method of claim 12, further comprising, based on 
the first, second, third, and fourth voltages, calculating an 
orientation of the receiver or of the transmitter. 

14. The method of claim 12, further comprising, based on 
the first, second, third, and fourth voltages, calculating an 
azimuthal angle of the receiver or of the transmitter. 

15. The method of claim 13, wherein calculating a dis- 
tance between the receiver and the transmitter and calculat- 
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ing an orientation of the receiver or of the transmitter are 
based on mutually decoupling the distance and the orienta- 
tion. 

16. The method of claim 15, wherein providing a receiver 
comprises only one receiver system. 

17. The method of claim 16, wherein calculating a dis- 
tance between the receiver and the transmitter and calculat- 
ing an orientation of the receiver or of the transmitter 
comprise receiving a location of either the receiver or the 
transmitter. 

18. The method of claim 12, wherein independently 
measuring each of the first, second, third, and fourth volt- 
ages is by receiving the quasistatic magnetic fields at a 
different frequency for the mutually orthogonal directions. 

19. The method of claim 12, wherein independently 
measuring each of the first, second, third, and fourth volt- 
ages is by receiving the quasistatic magnetic fields at a 
different time for the mutually orthogonal directions. 

20. The method of claim 13, wherein calculating a dis- 
tance between the receiver and the transmitter and calculat- 
ing an orientation of the receiver or of the transmitter are 
based on decomposing a near-ground magnetoquasistatic 
dipole into horizontal and vertical components. 

21. A method comprising: 

providing a receiver, the receiver comprising a first, 

second and third receiver coils configured to receive 
quasistatic magnetic fields in a first, second and third 
mutually orthogonal directions transmitted by a trans- 
mitter, the transmitter comprising a first, second and 
third transmitter coils and configured to transmit the 
quasistatic magnetic fields in a fourth, fifth and sixth 
mutually orthogonal directions; 

receiving the quasistatic magnetic fields through the 

receiver; 

measuring, with the receiver, a first voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the first receiver 
coil; 

measuring, with the receiver, a second voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a third voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the third receiver 
coil; 

measuring, with the receiver, a fourth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the first 
receiver coil; 

measuring, with the receiver, a fifth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a sixth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the third 
receiver coil; 

measuring, with the receiver, a seventh voltage originat- 

ing from the quasistatic magnetic fields due to coupling 
between the third transmitter coil and the first receiver 
coil; 

measuring, with the receiver, an eighth voltage originat- 

ing from the quasistatic magnetic fields due to coupling 
between the third transmitter coil and the second 
receiver coil; 
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measuring, with the receiver, a ninth voltage originating 
from the quasistatic magnetic fields due to coupling 
between the third transmitter coil and the third receiver 
coil; and 

based on the first, second, third, fourth, fifth, sixth, 

seventh, eighth and ninth voltages, calculating a dis- 
tance between the receiver and the transmitter, wherein 
the measuring of each of the first, second, third, fourth, 
fifth, sixth, seventh, eighth and ninth voltages is inde- 
pendent from the remaining voltages. 

22. The method of claim 21, further comprising, based on 
the first, second, third, fourth, fifth, sixth, seventh, eighth 
and ninth voltages, calculating an orientation of the receiver 
or of the transmitter. 

23. The method of claim 21, further comprising, based on 
the first, second, third, fourth, fifth, sixth, seventh, eighth 
and ninth voltages, calculating an azimuthal angle of the 
receiver or of the transmitter. 

24. The method of claim 22, wherein calculating a dis- 
tance between the receiver and the transmitter and calculat- 
ing an orientation of the receiver or of the transmitter are 
based on mutually decoupling the distance and the orienta- 
tion. 

25. The method of claim 24, wherein independently 
measuring each of the first, second, third, fourth, fifth, sixth, 
seventh, eighth and ninth voltages is by receiving the 
quasistatic magnetic fields at a different frequency for the 
mutually orthogonal directions. 

26. A method comprising: 

providing a receiver, the receiver comprising a first, 

second and third receiver coils configured to receive 
quasistatic magnetic fields in a first, second and third 
mutually orthogonal directions; 

providing a transmitter, the transmitter comprising a first, 

second and third transmitter coils and configured to 
transmit the quasistatic magnetic fields in a fourth, fifth 
and sixth mutually orthogonal directions; 

receiving the quasistatic magnetic fields through the 

receiver; 

measuring, with the receiver, a first voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the first receiver 
coil; 

measuring, with the receiver, a second voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a third voltage originating 

from the quasistatic magnetic fields due to coupling 
between the first transmitter coil and the third receiver 
coil; 
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measuring, with the receiver, a fourth voltage originating 
from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the first 
receiver coil; 

measuring, with the receiver, a fifth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a sixth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the second transmitter coil and the third 
receiver coil; 

measuring, with the receiver, a seventh voltage originat- 

ing from the quasistatic magnetic fields due to coupling 
between the third transmitter coil and the first receiver 
coil; 

measuring, with the receiver, an eighth voltage originat- 

ing from the quasistatic magnetic fields due to coupling 
between the third transmitter coil and the second 
receiver coil; 

measuring, with the receiver, a ninth voltage originating 

from the quasistatic magnetic fields due to coupling 
between the third transmitter coil and the third receiver 
coil; and 

based on the first, second, third, fourth, fifth, sixth, 

seventh, eighth and ninth voltages, calculating a dis- 
tance between the receiver and the transmitter, wherein 
the measuring of each of the first, second, third, fourth, 
fifth, sixth, seventh, eighth and ninth voltages is inde- 
pendent from the remaining voltages. 

27. The method of claim 26, further comprising, based on 
the first, second, third, fourth, fifth, sixth, seventh, eighth 
and ninth voltages, calculating an orientation of the receiver 
or of the transmitter. 

28. The method of claim 26, further comprising, based on 
the first, second, third, fourth, fifth, sixth, seventh, eighth 
and ninth voltages, calculating an azimuthal angle of the 
receiver or of the transmitter. 

29. The method of claim 27, wherein calculating a dis- 
tance between the receiver and the transmitter and calculat- 
ing an orientation of the receiver or of the transmitter are 
based on mutually decoupling the distance and the orienta- 
tion. 

30. The method of claim 29, wherein independently 
measuring each of the first, second, third, fourth, fifth, sixth, 
seventh, eighth and ninth voltages is by receiving the 
quasistatic magnetic fields at a different frequency for the 
mutually orthogonal directions. 
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